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TECHNICAL MEMORANDUM 


DESCRIPTION OF A MATHEMATICAL MODEL AND COMPUTER 
SIMULATION OF SEPARATION OF THE NOSE CAP FROM 
THE SOLID ROCKET BOOSTER 


I. INTRODUCTION 


The design of the decelerator subsystem for the Space Shuttle Solid Rocket Boosters includes a 
nose cap mounted at the forward end of the booster over the Drogue and Pilot parachute packs. The 
decelerator subsystem sequence of events begins with the ejection of the nose cap. In designing the sys- 
tem, it was necessary to determine what velocity increment had to be imparted to the nose cap so that it 
would separate without contacting the chute packs. Several different approaches were taken. A series of 
wind tunnel tests was conducted and a full scale test, using a rocket sled to propel a nose cap to high 
speed conditions, were run. In order to apply an analytic approach, a mathematical model was formu- 
lated. The equations of motion of the nose cap were derived by Steve Winder and Hughlen Murphree, 
both of the Systems Dynamics Laboratory, in 1974. The original computer program to implement the 
generation of nose cap trajectories was done by Hughlen Murphree. This report documents the model 
and computer program giving the system of equations, coordinate systems, and details of input and out- 
put data required to simulate a nose cap trajectory. 

The aerodynamic characteristics of the nose cap depend on whether or not it is in close proximity 
to the booster. When the cap is close to the booster the aerodynamic characteristics are functions of 
both angle of attack and displacement from the booster. When the cap moves wirhout a large afterbody, 
as it did in the rocket sled tests, the aerodynamic characteristics are functions of only angle of attack. 
(Mach number is essentially constant during the separation movement.) The program provides an option 
for simulation of either case. 

Because of design features of the SRB and the nose cap, the internal pressure in the cap could 
be very low at the time of separation. Thus, the external pressure would act to retard the separation. 
Provisions are made to include this effect by incorporating the internal pressu during separation as a 
function of separation velocity and displacement. 

The model incorporates two dimensions of translation ant! one of rotation. The British system of 
units is used throughout. 


II. THE MATHEMATICAL MODEL 


A. Coordinate Systems 

The basic coordinate system used is an inertial system with origin fixed at the center of gravity of 
the nose cap at time zero. The Xjj axis is horizontal and the Zjj axis is vertical with positive direction 

downward. Flight path angle, 7, is measured positive upward from the Xjj axis to the velocity vector. 
Attitude angle, 6 , is measured positive upward from the Xjj axis to the long axis of the nose cap, and 
angle of attack, a, is measured from the velocity vector to the long axis of the nose cap. The velocity 



components of the cap in this system, X^j and Zjj are called Ujj and W^, respectively. The velocity 
components of the booster are also given in this system and are called and Wjg. 

A body fixed coordinate system, with the origin fixed at the center of gravity of the nose cap, 
the Xj axis fixed on the long axis of the cap with positive direction toward the nose, and Zj axis per- 
pendicular to Xj, is used as the system in which the cap motion is described. The velocity in this system 
is denoted by Uj and Wj. These systems are illustrated in Figure 1. 

Finally, a coordinate system whose origin moves with the booster and whose axes are fixed 
parallel to X^j and Zjj is used to record the motion of points on the nose cap relative to the booster. 

These coordinates are denoted by X and Z with subscripts to identify the point being tracked. The 
velocity components of the cap in this system are called U and W. 


B. Equations of Motion and Integration 

A flat earth with constant gravity magnitude and direction is assumed. The booster motion 
during the time of interest is assumed to be of constant velocity magnitude and direction and the 
atmospheric density is assumed constant. These assumptions are warranted by the fact that the nose 
cap separates from the booster and clears the drogue pack in less than 0. 1 sec. 

The velocity and attitude of the mated booster and nose cap are given as Ujj, Wjj, and 9 , The 

AV of nose cap separation is given as DV. The first step is to express the velocity of the cap in the body 
fixed system. This is done by the equations: 


Uj = U u cos 0-W n sin 6 + DV 
Wi = Uji sin 9 + Wn cos 9 . 

Next using X m to represent nose cap mass and Xj for moment of inertia, the time derivatives of 
Uj, Wj, and 9 are given as: 

p + p 

Ui = -W, 0 - G sin 6 + — — ORIGINAL PAGE IS 

x m OF POOR QUALITY 

. . f n 

Wj = Uj 6 + G cos 6 + 

x m 


-F N (CG - CP) + x ml 
8 = ■■■ ■ 


0 



where 


G = acceleration due to gravity ORIGINAL PAGE IS 

= aerodynamic axial force POOR QUALITY 

Fjyj = aerodynamic normal force 

Fp = force due to difference between ambient pressure and nose cap internal pressure. 

CG = cap center of gravity location 

CP = cap center of aerodynamic pressure location 

X ml = Pitch moment coefficient times reference diameter times pV- A/2 . 

» 

These equations are numerically integrated by the Runge-Kutta method yielding Uj, Wj, and 0. Uj and 
Wj are transformed to U and W as follows: 

Ufj = Uj cos 9 + Wj sin v 


Wi i = — U i sin 9 + Wj cos 9 


u = u n -u IB 

W = W n -W 1B . 

* 

Finally, U, W, and 6 are integrated to produce X, Z, .and 0. 


C. Aerodynamics 

The aerodynamics of the nose cap are given in Reference 1 as functions of angle of attack, a, 
with separation distance from the booster as a parameter. If motion of the cap from the vicinity of 
the booster is to be simulated, the aerodynamic coefficients are determined by linear interpolation among 
points in tables with both oc and separation distance as independent variables. If motion of the cap 
without an afterbody present, as with the rocket sled tests, is to be simulated, the coefficients are deter- 
mined as functions of a alone using the data for separation distance equal to one diameter or greater. 

The differences between internal and external pressure are given as functions of separation 
distance for various cap separation velocities in Reference 2. This pressure difference applied over the 
base area of the cap is introduced in the equations of motion of the cap as a retarding force along the 
Xj axis. 
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D, Data Form 


Most applications of the nose cap trajectory model are done in terms of the angle of attack, 
dynamic pressure, altitude, and flight path angle pf the mated booster and nose cap at the instant of 
separation. These parameters do not appear explicitly in the model and were not used in the original 
computer program. The latest version of the program, however, has these parameters for input and 
incorporates the simple transformations that were originally done by hand each time a simulation was 
set up. One exception to this is altitude. The atmospheric density, p , which appears in the aerodynamic 
force and moment, and the gravitational acceleration, G, are functions of the altitude at which separation 
takes place. These must be calculated or determined from tabular data beforehand and must be input 
to the program. 


III. USING THE COMPUTER PROGRAM 


A. General Informat’on 

The source program is written in Fortran ’IV and is stored and functions on the Honeywell Sigma 
V computer. It is called NOSECAP. The data file is called DCAP. DCAP and NOSECAP are listed in 
Tables 1 and 2. 


B. Input 

The computer program is set up specifically for the aerodynamic data as it is given in Reference 
3. This reference gives axial force coefficient, normal force coefficient, and pitching moment coefficient 
as functions of angle of attack, a, with the parameter, separation distance divided by length, AX/D. 
Thirteen values of a are chosen and listed as AL(I). Axial force coefficients are specified as CAT1 (I) 
through CAT4(I) where the 1 signifies the value at separation distance, AX/D, equals 1 and 4 signifies 
separation distance equals 0. XCMTl(I) through XCMT4(I) similarly represents pitch moment coefflcient 
for separation distance 1 through 0. CNT1(I) through CNT6(I) represents the normal force coefficient 
since the reference gives six curves for normal force coefficient instead of four. If the number of entries 
in these tables is changed, the 'dimension statements must be changed accordingly. 

If the data are to be used to represent separation from the booster, the statement in line 146 of 
the program is used as part of the comment started in line 145, that is the letter C remains in the first 
column. If the sled test is to be simulated, the letter C is removed from the beginning of line 146 and 
the program proceeds with the quantity XCHK = 1 . 

The above data are part of the program and can cnly be changed by modification of the program. 
The following data are located in the data file and are read by the program. 

The difference between internal and ambient pressure is given in Reference 2 as a function of 
separation distance with nose cap velocity as a parameter. If the pressure difference is to be taken into 
account, the appropriate curve should be chosen and the data loaded as separation distance equals A01 
and pressure differentia] equals DPR. This table is dimensioned for twelve points. If more or fewer 
points are to be loaded, the dimension statement should be changed accordingly. If pressure difference 
is to be neglected, DPR should be loaded as zero. 
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The next data loaded are G, the gravitational acceleration in ft/sec**; RHO, the atmospheric 
density in slugs/fP; and XLENG, the length of the nose cap in inches. 

The next line of data is CPP, the location of aerodynamic center of pressure in percent of body 
length, and REFD, the reference diameter in inches. 

The next line of data is XI, the moment of inertia in slug ft"; ALPH, the angle of attack of the 

booster and cap before separation in degrees; DYNPR, the dynamic pressure before separation in lb/ft"; 
GAMMA, the flight path angle before separation in degrees; TH1D, the initial rate of change of attitude 
in degrees/sec; and DV, the separation velocity increment in ft/sec. 

The next line of data is WT, the weight of the cap in lb; CGP, the location of the center of 
gravity of the cap in percent of body length; XCG, the X coordinate of the center of gravity of the cap, 
in body fixed coordinates (usually zero); ZCG, the Z coordinate of center of gravity of the cap in body 

fixed coordinates (usually zero); and AREA, the reference area of the cap in ft^. 

The last line of input is XCONE, the distance from the cap CG to tire base of the cap in feet; 
and ZCONE, the distance in feet from the center line of the nose cap to a point on the base of the cap. 
XCONE and ZCONE establish the points P and Q which are tracked in the trajectory printout (Fig. 1 ). 
This completes the data file. 

Besides printing the trajectory of nose cap CG, and points P and Q as functions of time, the 
computer program has a subroutine called DIST which tracks an extreme point on the nose cap (usually 
the comer of the thruster post) and determines the distance up the side of the drogue pack where this 
point intersects the drogue pack envelope (Fig. 2). This distance is printed between time steps where 
the crossing of the pack envelope occurs. This subroutine contains constants that could be changed if 
a different cap configuration is studied. These constants are illustrated in Figure 2. Note that the angle 
PP is used in the subroutine in radians, ’ 

The integration time step may be adjusted if a different cap configuration is studied (for instance 
if wind tunnel tests are simulated). The time step DELTA is specified on line 130 of the main routine. 

Also, the value of DELTA is-used in line 289. These values should always be equal. The value, 2 , has 
been found appropriate for simulating motion of the full scale SRB nose cap. For simulation of wind 
tunnel tests where the scale factor is 1/8, for example, the integration time step should be scaled by a 

factor %/ 1/8, 

The time point at which to stop the computation is given at line 87 as TSTOP. At present, the 
value TSTOP = 0.16 is used. 


C. Output 

The input data file is printed out in the first block of output. The printout then proceeds with 
trajectory parameters at regular time intervals. A list of the mnemonics and the parameters they repre- 
sent is given in Table 3. 
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IV. CONCLUSION 


Completion of the sled test presented an opportunity to compare the nose cap separation simula- 
tion with test results. Motion picture film of the test was analyzed and a graphic representation of the 
nose cap motion was drawn. A simulation was then made using the test conditions and a similar graph 
was overlaid on the test results. This is shown in Figure 3. The coordinates used are those of the film 
analysis. The abscissa is along the base of the nose cap before separation. The origin is at the leading 
edge of the nose cap before separation. Positions of the nose cap are shown for both test results (solid 
line) and simulation (broken line) at various times. The times do not coincide exactly, in general, 

because the simulation time step was 2’® sec and the film speed was 0.0125 sec per frame. Nevertheless, 
the simulated motion is very similar to the observed motion. 

The test condition was the design limit condition of 80 deg angle of attack and dynamic pressure 

200 lb/ft-. The data file shown in Table 1 is the one used for this simulation and all of the test condi- 
tions are included there. This particular simulation serves as a sample of the program function. It is 
identified as run 0464 and is kept in the flies of the Dynamics and Trajectory Analysis Branch. 

The good agreement between the test and simulated data has established the validity of use of the 
simulation to analyze the separation motion. 
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Figure 1. Coordinate systems and preload parameters. 
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Figure 3, Comparison of motion of base of nose cap as simulated and observed in sled test. 



TABLE 1. LISTING .OF DCAP 


* 


i 0,0,5. 0 , i0.0 # 

2.000 30 . 0 , 35 . 0 , 40 . 

3.000 0 ., 0 ., 0 ., 0,,0 
* 4.000 0 ., 0 ., 0 ., 0.,0 

5.000 32 . 3 , , 0 < y 2 , 75 . 

6.000 . 667 , 68.03 

7.000 52 . , 80 . , 200 , , 

8.000 300 ,, . 66 , 0 . ,0 

9.000 2 . 125 , 2.83333 


15. 0, f §0.-0, 25. 0, 
0/45.0,50,0,300.0 
* / 0 . , 

» / 0 » > 


0 • / 0 . , 76 ♦ , 
., 25 * 24 , 


P 


OP’GINAL PAGE IS 
Of POOR QUALITY 
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TABLE 2. LISTING OF NOSBCAP 


1.000 C 

2.000 

3.000 
L( 13 ) 

4 . 000 

5.000 

6 .000 

7 .000 

8 . 000 
9.000 

10.000 
11 . 000 
12.000 

13.000 

14.000 

15.000 
-1 . 3, 

16.000 

17 . 000 
-1 .3, 

18.000 

19 . 000 
,-1.3, 

20 .000 
21 . 000 

1 > “ 1 * 3 , 
22.000 

23.000 
• » . 04 , 

24 .000 

25.000 
-. 02 , 

26.000 

27.000 
3, 

28.000 

* 


NOSE CAP SEPARATION STUDIES 
REAL M, Ml 

D I HENS I ON CNTl(13),XCMTlCi3),CATl < 13 ) , CNTS ( 13 ) , CNT6 ( 13 ) , A . 

DIMENSION CNT2 ( 1 3 ),XC(1T2( 13 ), CAT 2(13 ) 

DIMENSION CNT3(13 ),XCMT3( 13 ),CAT3(13 ) 

DIMENSION CNT4( 13 ),XCMT4< 13 ) , CAT4 ( 13 ) 

DIMENSION AOl ( 12), DPRC 12 ) 

COMMON/SEC D/ T H i D D , DUM0<59) 

C0MM0N/FUNT/Y(3),DUM1 (57 ) 

COMMON/FIRT/F (3 ),DUM2(57) 

C0MM0N/LAB/U1P,TH1P,REFD 
COMMON XCONE , ZCONE 

EQUIUALENCE (TH1D,F<1)),(U1D,F<2)‘>,<U1D,F<3)> 

EQUIUALENCE <THi,Y(l)),<Ul,Y(2))(Ui,Y<3)) 

DATA (CATl(I), 1-1,13)/. 4, .42, .38, .25,.!,. 02, .22, -.18, -.8, 


*-l .38,-1. 38,-1 .4/ 

DATA ( CAT2 ( I ), 1-1,13)/. 4, .42, .3,. 18,0. ,-.06,. 16, -.19, -.9, 


*-l .38,-1 .38,-1 . 4/ 

DATA ( C AT3 ( I ), 1-1,13)/. 4,. 42,. 28,. 1,-. 02, -.16,. 04,-. 2, .-1. 


*-l .38,-1 . 38,-1 .4/ 

DATA (CAT4(I ), 1-1,13)/. 4, .42, .21, .06,-.!, -.18,0., -.21,-1. 


*-l .38,-1 .38,-1 .4/ 

DATA (XCMT1 ( I ), 1-1, 13)/0. . .02, .03, . 02 , - . 02, - . 05 , - . 1 , - . 1 , 0 
*.05, .04,0./ 

DATA ( XCMT2 ( I ),I«l,13)/0.,.02,.03,0.,-.03,-.07,-.13.-.12, 
*.04, .05, .04,0./ 

DATA (XCMT3(I ), I-l,13)/0., .01, . 02, -.02, -.04, -.09,-. 15, -.1 
*-.03, .03, .04, .04,0./ 


ORIGINAL PAGE 19 
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TABLE 2. (Continued) 


£9.000 DATA < XCMT4 < I ) , I ■ 1 , 1 3 )/0 . , . 0 1 , . 01 , - . 04 , - . 08 , - . 1 2 , - . 18 , - . 

14," .06, 

30.000 * .01 , .03, .04,0./ 

31 .000 DATA ( Al( I ), 1-1, 13 )/0. , 15. ,30 . ,45. ,60. ,75. ,90. , 105. , 120. , 
135. , 

32.000 *150. ,165. , 180./ 

33.000 DATA ( CNT 1 ( I ),I-l,13)/0., .38, .8, *82, .8, .7, .55, .44, .33, .12 


34.000 

35.000 

13, 

36.000 

37.000 

14, 

38.000 

39.000 
.15, 

40.000 

41 .000 

37. . 16, 

42.000 

43.000 

39. . 17, 

44.000 

45.000 

46.000 

47 . 000 

48.000 

49.000 

50.000 

51.000 

52.000 

53.000 

54.000 

55.000 

56 . 000 

* 


*-.08,-. 1,0. > 

DATA ( CNT2 ( I ) , I » 1 , 1 3 ) /0 . , . 39 , . 82, . 9, . 86 , . 76 , . 57 , . 45 , . 34 , . 
*- . 02 , -. 1 , 0 ./ 

DATA (CNT3(I),I-l,13)/0., .4, .84, .98, .96, .81, .61, .46, .35, • 

a- . 02 , - . 1 , 0 ./ 

DATA ( CNT4 < I ), 1-1,13)/®., .41, .9, 1.01, 1.02, .9, .65, .49, .36, 
*- . 02, - . 1 , 0 ./ 

DATA ( CNT5(I ) , I - 1 , 13 ) /0 . , . 42, . 97, 1 . 1 3, 1 . 1 6, 1 . 05 , .72, .51, . 
*- . 02, - .1,0./ 

DATA ( CNT6 ( I ) , I - 1 , 13 )/0 . , . 43 , 1 . 12 , 1 . 35 , 1 . 46, 1 . 4 , . 8? , . 53, . 

*-.02, -.1,0./ 

25 READ (5,92 ), (AOl ) 

READ (5,92), (DPR) 

URITE (6, 92), (AOl) 

URITE (6, 92), (DPR) 

READ (5,2), (G,RHO,XLENG) 

READ (5,2), (CPP,REFD) 

READ (5,92 ), ( XI , ALPH, DYNPR, GAMMA, TH1D, DU ) 

READ (5,2), ( (JT, CGP, XI , Z1 , AREA ) 

READ (5,2), ( XCONE , ZCONE ) 

XLEN-XLENG 

ANGLEA - ATAN(XLENG/(REFD*.5 ) ) 

SIDEL - XLENG/ N ( ANGLEA ) 
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TABLE 2. (Continued) 


57*000 SIDEL-SIDEl/12. 

58.000 U1P » DU 

59.000 TH1P - ALPH+GAMMA 

60.000 URITE(6,92), (G,RHO,XLENG) 

61.000 URITE(6,92 ), (CPP,REFD ) 

62.000 URITE(6,92) , ( XI , ALPH, DYNPR, GAMMA, TH1D, DU ) 

63.000 UR I TE ( 6 , 92 ) , < UT, CGP . XI , Z 1 , AREA ) 

64.000 URITE(6,92)» (XC0NE,ZC0NE) 

65.000 92 FORMAT ( 6E15.6) 

66.000 C-57. 29578 

67.000 ALPH-ALPH/C 

68.000 GAMMA-GAMMA/C 

69.000 TH1 -ALPH+GAMMA 

70.000 UUU-SGRT(2.0*DYNPR/RHO ) 

71.000 U1 I -UUU*COS(GAMMA ) 

72.000 UlI--UUU*SIN(GAf1MA ) 

73.000 U1B-U1I 

74.000 U1B-U1I 

75.000 TH1DD-0. 

76.000 PSI-ATAN2(ZC0NE,XC0NE) 

77.000 TH1D - TH 1 D/C 

78.000 TH11-TH1 

79.000 THETA - TH1 

80.000 U1 - U1I*C0S(THI) - Ul I #SIN ( TH1 ) + DU 

81.000 U1 - U1I*SIN(TH1> + U 1 I *COS ( TH1 ) 

82.000 C THE ABOUE EQUATIONS PUT U1 AND U1 IN THE BODY AXIS SYSTEM 

83.000 III-0 

84.000 K • 0 

85.000 L ■ 1 2 

86.000 N-13 

87.000 TSTOP-0.16 

88.000 XM-UT/G 

89.000 THA-PSI-TH1 

90.000 R-8. *0.46/12. 
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TABLE 2. (Continued) 


91.000 
9B.000 

93.000 
9-4.000 

95.000 

96.000 

97.000 

98.000 

99 . 000 
100.000 
101 . 000 

132.000 

103 . 000 

104 . 000 

105.000 

106.000 

107.000 

108.000 

109.000 

110. 000 
111.000 
112.000 

113. 000 

114. 000 

115. 000 

116. 000 

117.000 

118.000 

119.000 

120 . 000 
121 .000 
122.000 

123.000 

124.000 


TLEN- (2.56*8. )/12. 

FLEN« (8. 1*8 )/12. 

THI--TH1 

XA«X1~XC0NE*C0S(THI )-ZCONE*SIN(THI )+R*SIN(THI ) 
ZA-Z1-XC0NE*SIN(THI ) + ZCONE*COS(THI )-R*COS ( THI ) 
G1 -108./57.3-THI 
G2- (72.Z57.3 )+THl-(90 . /57.3) 

M “ TAN ( G2 ) ■■ 

B»ZA-n*XA 

XB-X1-XC0NE*C0S(THI ) +ZCONE*S IN ( TH I )-R*SIN(THI ) 
ZB-Z1-XC0NE*SIN(THI )-ZCONE*COS < THI )+R*COS(THI ) 
ACl-fUM+1 . 

BC1 ■-2.*XA-2.*t1*ZA + 2.*B*M 

CC1»XA*XA+ZA*ZA-2.*B*ZA+B*B-FLEN*FLEN 

XA1- (-BCi + SQRT(BCl*BCl-4 .*ACl*CCi ) )/(2.*ACl ) 

Z A 1 * <ri*XAi+B) 

Ml- (ZB-ZA )/(XB-XA) 

B2-ZA1-H1TXA1 
AC2«imm + i . 

BC2--2.£XA1+2.*J11*B2-2.*I11*ZA1 

CC2-XA1*XA1+B2*B2-2.*ZA1*B2+ZA1*ZA1-TLEN*TLEN 
XBlr <-BC2+SQRT(BC2*BC2-4. *AC2*CC2 ) )/<2.3AC2 ) 

ZB 1 ■ <M1*XB1 + B2) 

PI-TAN ( G1 ) 

XTP-XA 

ZTP-ZA 

URITE(3 ),XTP,ZTP 

XTP-XB 

ZTP m ZB 

URITEC3 ),XTP,ZTP 

XTP-XB1 

ZTP-ZB1 

URITE(3),XTP # ZTP 

XTP-XA1 
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TABLE 2. (Continued) 


125.000 

126.000 

127.000 

128. 000 

129.000 

130.000 

131.000 

132.000 
i 33.000 

134 . 000 

135.000 2 

136.000 200 

137.000 500 

138.000 

139.000 

140.000 501 

141.000 502 

142 . 000 

143.000 C 

144.000 C 

145.000 0 

146.000 C 

147.000 C 

148.000 C 

149.000 

150.000 C 

151.000 901 

152.000 

153.000 

154.000 

155.000 C 

156.000 902 

157.000 

158.000 


ZTP-ZA1 

URITE(3),XTP f 2TP 

XTP-XA 

ZTP-2A 

UR ITE ( 3 ),XTP,ZTP 
DELT- .00390625 
TIME-0. 0 
PRTIM - 0. ; 

PLTSAU - 0. 

I - 1 

FORMAT (5E15.6) 

IND-1 

CONTINUE 

ALPHA-XTAN2(U1,U1 ) 

IF (ALPHA ) 50 1 , 502 * 502 
ALPHA-360. 0+ALPHA 
CONTINUE 

XCHK-X1/ (REFD/12. ) 


THIS IS ONLY USED FOR THE SLED TEST EUALUATIONS 
XCHK-1 .0 


IFCXCHK.GE.l. ) GO TO 900 

THIS IS FOR XCHK-0.3 

CALL TBL(ALj,CNT2,ALPHA # N,CN2) 

XCHK1 -0 . 3 

CNC1-CN2 

IF (XCHK.GE.0.3 ) GO TO 900 
THIS IS FOR XCHK-0.2 
CALL TBL(AL,CNT3 # ALPHA,N,CN3) 
IF(XCHK.GE.0.2) GO TO 802 
XCHK1-0.2 



TABLE 2. (Continued) 


o\ 


159 .000 

160.000 
161.000 
166.000 

163.000 
16-1.000 C 

165.000 

166 . 000 

167 . 000 

168.000 

169. 000 

170.000 

171.000 

172.000 

173.000 C 

174.000 

175.000 

176.000 

177.000 

178.000 

179.000 

180 . 000 
181 . 000 

182.000 C 

183.000 

184 . 000 

185 . 000 

186.000 

187.000 

188 . 000 

189.000 

190.000 C 

191 . 000 

192.000 

* 


CNC1-CN3 
GO TO 903 

802 XCHK2-0.2 
CNC2-CN3 
GO TO 909 

THIS IS FOR XCHK-0.1 

903 CALL TBL(AL,CNT4,ALPHA,N,CN4) 
I F ( XCHK . GE . 0 . 1 ) GO TO 803 
XCHK1-0. 1 

CNC1 »CN4 
GO TO 904 

803 XCHK2-0.1 
CNC2-CN4 
GO TO 909 

THIS IS FOR XCHK *0.05 

904 CALL TBL(AL,CNT5, ALPHA, N,CN5 ) 
IF (XCHK. GE. 0*05) GO TO 804 
XCHK 1*0. 05 

CNC1 -CN5 
GO TO 905 

804 XCHK2 ■ 0.05 
CNC2-CH5 
GO TO 909 

THIS IS FOR XCHK • 0 . 0 

905 CALL TBL(AL,CNT6, ALPHA, N,CN6 ) 
IF (XCHK. GE. 0.0) GO TO 805 
CN-CN6 

GO TO 911 

805 XCHK2-0.0 
CNC2-CN6 
GO TO 909 


THIS IS FOR XCHK *1.0 
900 CALL TBL< AL.CNT1 , ALPHA, N,CN1 ) 
XCHK 2 • 1 . 0 




TABLE 2. (Continued) 


193.000 



CNC2-CN1 

194.000 



IF ( XCHK . GE . 1 • 0 ) GO TO 910 

195.000 


909 

CN-( (XCHK-XCHIC1 )/< XCHK1-XCHK2 ) )* ( CNC 1 -CNC2 ) +CNC1 

196.000 



GO TO 911 

197.000 


910 

CN-CN1 

198 . 000 


911 

CONTINUE 

199.000 



IF ( XCHK . GE * 1 . ) GO TO 700 

200.000 

C 


THIS IS F OR XCHK *0-2 

201 .000 


701 

CALL TBLCAL, XCMT2, ALPHA, N,XCM2 ) 

202 . 000 



CALL TBL<AL,CAT2, ALPHA, N,CA2 ) 

203 . 000 



XCHK1-0.2 

204 . 000 



XCtll-XCME 

2O5.000 



CA1-CA2 

206.000 



IF ( XCHK • GE . 0 . 2 ) GO TO 700 

207.000 

C 


THIS IS FOR XCHIC-0.0S 

208.000 


702 

CALL TBL ( AL , XCNT3, ALPHA, N, XCP13 ) 

209.000 



CALL TBL<AL,CAT2, ALPHA, N,CA3 ) 

210.000 



IF ( XCHK . GE .0.05) GO TO 601 

211 . 000 



XCHK 1*0. 05 

212.O00 



XCF11-XCH3 

213.000 



CA1-CA3 

214.000 



GO TO 703 

215 . 000 


601 

XCHK2-0 . 05 

216.000 



XCM2-XCM3 

217.000 



CA2-CA3 

218.000 



GO TO 709 

219.000 

c 


THIS IS FOR XCHK-0.0 

220.000 


703 

CALL TBL(AL, XC11T4, ALPHA, N, XCM4 ) 

221 . 000 



CALL TBL(AL,CAT4, ALPHA, N,CA4 ) 

222.000 



IF(XCHK.GE.0.0) GO TO 605 

223.000 



XCN-XCM4 

224.000 



CA-CA4 

225.000 



GO TO 711 

226.000 


605 

XCHK2-0.0 



TABLE 2. (Continued) 


227.000 


XCM2-XCM4 

228 .000 


CA2*CA4 

229.000 


GO TO 709 

230.000 

C 

THIS IS FOR XCHK-1.0 

231.000 

700 

CALL TBUAL, XCNT1, ALPHA, N,XCM6 ) 

232.000 


CALL TBL(AL,CAT1, ALPHA, N,CA6 ) 

233.000 


XCHK2-1 .0 

234 .000 


XCN2-XCM6 

235.000 


CA2-CA6 

236.000 


IFCXCHK.GE.l. ) GO TO 710 

237.000 

709 

XCP1»( ( XCHK-XCHK1 ) / ( XCHK1-XCHK2 ) )* ( XCH1 -XCN2 ) + XCPU 

238.000 


CA- ( (XCHK-XCHKi )/ ( XCHK1 -XCHK2 ) )* ( CA1-CA2 ) + CAl 

239.000 


GO TO 711 

240.000 

710 

XCM-XCM6 

241.000 


CA-CA6 

242.000 

711 

CONTINUE 

243.000 


AA- 0,S*RHO*(U1*U1 + U1*U1 )*AREA 

244.000 


XMl«AA*XCn*(REFD/12.0 > 

245.000 


C1P"-(XLEN/12.0 )X(CGP-CPP ) 

246.000 


FA --AAXCA 

247.000 


FN - -AAXCN 

248.000 


XMOMEN-TH1DDXXI 

249.000 


TH1DD- ( ( C1PXFN )+XHl )/XI 

250.000 


U ID ■ ( 1 . 0/XM)*(XM*Ul*THlD + XN*GXCOS(THl )+FN) 

251 . 000 


U1D- ( 1 .0/XM )X(-XMXUl*THlD-Xn*G*SIN(THl ) + FA + FP ) 

252.000 


IF C TINE-0. 0 ) 600, 600, 41 

253.060 

41 

CONTINUE 

254.000 


CALL RUNGEtlND, TIME, DELT, 3,1 ) 

255 . 000 


IND-IND+1 

256.000 


IF( IND-4 >500,500,600 

257.000 

600 

CONTINUE 

258.000 


U1I - U 1 XC05 ( TH 1 ) +U1*SIN(TH1) 

259.000 


U1I - -U1XSINCTH1) + U1XCOS (TH1 ) 

260.000 


FP A - ATAN2(U1I,U1I )XC 



TABLE 2 . (Continued) 


261.000 

262.000 

263.000 

264 . 000 


U 1 1 • U1I - U 1 B 
U1I • U1I - U1B 
UEL-SQRT(U1I*U1I + 1)1I*U1I ) 
URITE(6,39 ) , T IME 

t 

URITE(6,3 ), (ALPHA ) 

t 

URITE(6, 1001 ) , <FN,FA i 

265.000 


266.000 


267.000 


URITE(6,4 )/<XM0f1EN) 
♦ 

268.000 


* 

WRITE (6, 1002 ), (U1I,U1I,FC1) >,FPA 

t 

URITE(6,2222), (UEL) 

♦ 

269.000 


270.000 


l 

URITE ( 6, 8 ) , XCM 

271.000 

8 

4 

FORMAT C1X,4HXCM-,E15.6) 

t 

URITE<6,7), AA, XM1 ,C1P,U1,U1,CN 

272.000 


273.000 

7 

• 

FORMAT (IX, 4H A A- , E 1 5 . 6, 4HXM1 - , E 15 . 6, 4HC IP- , E15 . 6, 4H 

E15.6, 


i 

274.000 


*4H U1-,E15*6,4H CN-,E15.6) 

f 

275. 00‘0 


4 

TH1DS-F ( 1 ) 

276.000 


TH1S-VC 1 ) 

277.000 


TH1S1-TH1S*C 

278.000 


U1S- V(2 ) 

279.000 


U1S-V(3 ) 

280.000 


F ( 1 ) -0 . 0 

281 . 000 


F ( 2 ) -U1 1 

282.000 


F ( 3 ) • U1 I 

283.000 


V( 1 ) -0.0 

284.000 


Y ( 2 ) - 2 1 

285.000 


V ( 3 ) - X 1 


ORIGINAL PAGE IS 
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u 

o 


TABLE 2. (Continued) 


286.000 

287.000 
288.060 

289.000 
290.0C0 

291.000 

292.000 

293.000 

294.000 

295.000 

296.000 

297.000 

298.000 

299.000 

300 . 000 

301.000 

302.000 

303 . 000 

304.000 

305.000 

306.000 

307.000 

308.000 

309.000 

310 . 000 

311 .000 

312.000 

313.000 

314.000 

315 . 000 

316.000 

317.000 

318.000 

319.000 

* 


IND-1 

IF(TIME-0,0 >800,800,300 
300 CONTINUE 

TINE-TIME-. 00390625 

400 CONTINUE 

CALL RUNGE< IND, TIME, DELT, 3,0 ) 
IND-IND+1 

IF ( IND-4 1400,400,800 
800 TH 1 -THIS 
2 1 * Y ( 2 ) 

XI - Y<3 ) 

XI I- (X1XC0SCTH11 )-Zl*SIN(THU ) )*12. 
CALL TBL(A01,DPR,X11,L,DPRE) 
FP--DPRE*AREA*144. 

URITE(6, 1003 ) , (Y(3),Y(2),TH1S1,FP) 

THA -PSI-TH1 

TH3-TH1 

XP - XI - XC0NE*C0S ( -TH3 ) 

2P - 22 - XC0NE*SIN(-TH3 ) 

XQ - XI - XC0NE*C0S(~TH3 > 

ZQ - Z1 - XC0NE*SIN(-TH3 ) 

IF (TIME-. 08 1401,402,402 

401 CONTINUE 
XTP-XQ 
ZTP-ZQ 

IJRITEd ) , XTP , ZTP 
XTP-XP 
ZTP-ZP 

URITEt 1 ), XTP, ZTP 

402 CONTINUE 
URITE(6,6 ), ( XO , ZQ ) 

C 

UR I TE (6,5) , (XP, ZP 1 
URITE(6, 1013 ), (JIB, DELT 
F ( 1 ) ■ TH1DS 


ZC0NE*SIN(-TH3> 
ZC0NE*C0S(-TH3 ) 
ZCONE* SIN ( -TH3 ) 
ZC0NE*C0S(-TH3 ) 


TJ 

l 



TABLE 2. (Continued) 


320.000 

321 .000 

322.000 

323.000 

324.000 

325.000 

326.000 

327.000 

328.000 

329.000 

330.000 

331.000 

332.000 

333.000 

334 . 000 

335.000 

336.000 

337.000 

338.000 

339.000 
.6. 10X> 

34e.000 

341.000 
15.6, 

342*000 

343.000 

344.000 

345.000 

346.000 

347.000 

348.000 

349.000 

350.000 

* 


V ( 1 ) »TH1S 
Y(2 >«U1S 
Y ( 3 ) - U1 S 
I - I + 1 
490 CONTINUE 

IF (TIME .LE. 0.0195313) GO TO 201 
CALL DIST(TH11,TH1S,X1,Z1,K) 

201 CONTINUE 
XS1-X1 
ZS1-Z1 
THS-TH1 

IF<TIME-0.0)41,41,49 
49 CONTINUE 

I F < T I ME-TSTOP >24,251,251 
24 CONTINUE 
GO TO 200 
251 CONTINUE 
GO TO 25 

1001 FORMAT < 1H,2HFN,3X,E15.6,2X,2HFA,3X,E15.6 ) 

1002 FORMAT ( 1H, 2HU , 3X, El 5 * 6, 2X, 2HU , 3X, E15 . 6 , 2X, 3HTHD, 3X, El 5 


QO 
-n x 

x a 
O 5 

o > 
» F 

o x 
> § 

I- 

c n 


*20HFLIGHT PATH ANGLE - ,E15.6) 

1003 FORMAT ( 1H, 3HXCG, 2X, El 5 . 6 , 2X, 3HZCG , 2X, E15 . 6, 2X , 3HTH 1 , 3X, E 


*2X,2HFP,3X,E15.6) 

3 FORMAT ( 1H,5HALPHA,3X,E15.6 ) 

4 FORMAT (1H,6HM0MENT,3X,E15.6) 

5 FORMAT ( 1H,2HXP,3X,E15.6,2X,2HZP,3X,E15.6 ) 

6 FORMAT C1H,2HXG,3X,E1S.6,2X,2HZQ,3X,E15.6 ) 
2222 FORMAT ( 1H, 3HUEL, 3X, E15 . 6 ) 

39 FORMAT ( 1H, //IX, 4HTIME,3X, E15.6 ) 

1013 FORMAT ( IX, 5HU1B" , E15 . 8, 10X, 5HDELT- , E15 . 8 ) 
END 
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TABLE 2. (Continued) 


351.000 

352.000 

353.000 

354.000 

355.000 

356.000 

357.000 
0 ), 

358.000 

359 . 000 

360.000 

361.000 

362.000 

363 .000 

364 . 000 

365 .000 

366.000 

367.000 

368.000 

369.000 

370 . 000 

371 . 000 

372 .000 

373.000 

374.000 

375.000 

376.000 

377.000 

378.000 

379 . 000 

380.000 

381 . 000 

382.000 

383.000 

384 . 000 

* 


SUBROUTINE RUNGE ( KUTT A, TI (IE , DT , NUAR, NDUAR ) 

C KUTTA IS A CONTROL INTEGER 
C TINE IS TINE OF INTEGRATION 
C DT IS THE TINE INCREMENT 

C NUAR - NUMBER OF UARIABLES TO BE INTEGRATED ONCE 
C NDUAR- NUMBER OF UARIABLES TO BE INTEGRATED TUICE 

DIMENSION Cl( 60 ) , C2 ( 60),C3< 60),C4( 60),CD1( 60),CD2< 6 

1 CD3 ( 60 ) , CD4 ( 60 ) , SX < 60),SXD< 60) 


1 

COMMON/FUNT/ X< 60) 
COMMON/FIRT/ XDC 60) 
COMMON/SECD/XDD(30 ) 

GO TO (1,2, 3, 4), KUTTA 
DO 10 I - 1 , NUAR 



SX(I ) - X ( I ) 

Cl ( I ) - XD( I )*DT 

3l 

10 

X(I) - SX ( I ) + 0 . 5SC 1 ( I ) 


IFCNDUAR.EQ.0 )GO TO 50 


DO 100 1-1, NDUAR 

S? 


SXD(I) - XD ( I ) 

iC "0 


CD 1 ( I ) - XDD(I)*DT 

c 5 
0 

100 

XD ( I ) - SXD(I) + 0.S*CDKI) 

r pi 

50 

CONTINUE 



TIME - TIME + 0.5*DT 
RETURN 

2 DO 20 I - 1 , NU AR 
C2 ( I ) - XD(I)*DT 

20 X(I) - SX < I ) +0 . 5*C2 ( I ) 
IFCNDUAR.EQ.0 )GO TO 51 
DO 200 I - 1 , NDUAR 
CD2CI) - XDD ( I )*DT 
200 XD ( I ) - SXD< I ) + 0.5*CD2(I ) 
51 CONTINUE 
RETURN 

3 DO 30 I • 1 , NUAR 



TABLE 2. (Continued) 


385. 0O0 


C3 ( I ) ■ XD(I)*DT 

386.000 

30 

X ( I ) • SX ( I ) + C3 ( I ) 

387.000 


IF<NDUAR. EQ.0 )G0 TO 52 

388.000 


DO 300 I ■ 1 > NDUAR 

389.000 


CD3 ( I ) » XDD(I)*DT 

390.000 

300 

XD ( I ) » SXD ( I ) + CD3 ( I ) 

391 . 000 

52 

CONTINUE 

392.000 


TIME - TIME±0.5*DT 

393.000 


RETURN 

394.000 

4 

DO 40 I-i,NUAR 

395.000 


C4 ( I ) - XD(I)*DT 

396.000 

40 

X(I) - SX(I)+(Cl(I)+C4(I)+S.*(C2(I)+C3CI)))/6.0 

397.000 


IFtNDUAR .EQ.0 )GO TO 53 

398.000 


DO 400 I - 1 , NDUAR 

399.000 


C D4 ( I ) - XDD(I)*DT 

400 . 000 

400 

XD(I) * SXDC I )+(CDl ( I )+CD4( I )+2.0*(CD2(I ) + CD3< I ) ) )/6.0 

401.000 

53 

CONTINUE 

402.000 


RETURN 

403.000 


END 


to 

CO 
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TABLE 2 . (Continued) 


404,000 


SUBROUTINE TBL ( X , V , XI , N , Vi ) 

405.000 C 


THE DATA MUST BE IN ASCENDING 

406 . 000 


DIMENSION X(i ), YC1 ) 

407.000 


DO 93 1*1, N 

408.000 


I F ( X ( I )-Xl >93,95,91 

409.000 

91 

DX • XI -X ( 1-1 > 

410.000 


DX1 *X( I )-XCI-l ) 

411. 000 


D V • V ( I > - V C 1-1 ) 

412.00O 


V1-V(I-1 ) + D **(DX/DX1 ) 

413.000 


GO TO 94 

414.000 

95 

VI «V( I ) 

415,000 


GO TO 94 

416.000 

93 

CONTINUE 

417.000 

94 

CONTINUE 

418.000 


RETURN 

419.000 


END 


ORDER 



TABLE 2. (Continued) 


420.000 


FUNCTION XTAN2( A, B ) 

421.000 


IF ( B .EQ. 0.0) GO TO 

422.000 


XTAN2» ATAN2 ( A, B )*57.3 

423.000 


RETURN 

424.000 

100 

IF( A >10, 16, 15 

425.000 

10 

XT AN2“ -90 • 0 

426.000 


RETURN 

427.000 

16 

XTAN2-0.0 

428.000 


RETURN 

429.000 

15 

XT AN2 * 90 . 0 

430.000 


RETURN 

431.000 


END 


to 

Ol 


$3 

S| 

§ P 





t J 
CK 


TABLE 2. (Concluded) 


432.000 


SUBROUTINE DIST(THll,THi$,Xl,Zl,K) 

433.000 


IFCK .EG. 1 ) GO TO 330 

434.000 


K-0 

435.000 


XX 1 ■ -2 . 552083 

436.000 


ZZ1 -2.5&S367 

437.000 


XX2--2. 125 

438.000 


ZZ2-2. 516667 

439 . 000 


PP- .3131086 

440 .000 


S0-SIN< TH11 ) 

441 . 000 


CO-COS ( TH1 1 ) 

442.000 


XX3-XX2*C0+ZZ2*S0 

443.000 


ZZ3«ZZ2*C0-XX2*S0 

444 . 000 


S - S I N ( TH 1 S ) 

445 . 000 


C«C0S(TH1S ) 

446.000 


XXS-X1+XX1*C+ZZ1*S 

447 . 000 


ZZ5-Z1+ZZ1*C-XX1*S 

448 . 000 


XX6 - XX5-XX3 

449 . 000 


ZZ6 * ZZ5-ZZ3 

450.000 


PP0--TH11-PP 

451.000 


S 1 - 5 IN ( PP0 ) 

452.000 


Cl -COS ( PP0 ) 

453.000 


XX?-XX6*C1+ZZ6*S1 

454 . 000 


ZZ?-ZZ6*C1-XX6*S1 

455.000 


IF ( ZZ7 .LT. 0. ) GO TO 340 

456 . 000 


XX8-XX7 

457.000 


ZZ8-ZZ7 

458.000 

330 

RETURN 

459.000 

340 

XXX- (XX7-( (ZZ7*<XX8-XX7 ) )/<ZZ8-ZZ7)) )*12. 

460.000 


UR I TE ( 6 , 360 ) » (XXX) 

461 .000 

360 

FORMAT ( ///*7H DISTIN, 3X,E15.6,/// ) 

462.000 


K- 1 

463.000 


RETURN 

464.000 


END 


«§ 

31 

S? 

g? 

> s 

52 



TABLE 3. OUTPUT PARAMETERS 


Mnemonic 

Definition and Units 

TIME 

Time (sec) 

ALPHA 

Angle of attack (deg) 

FN 

Normal force (lb) 

FA 

Axial force (lb) 

MOMENT 

Product of angular acceleration and moment of inertia (lb/ft) 

U 

Velocity relative to booster component in X direction (ft/sec) 

W 

Velocity relative to booster component in Z direction (ft/sec) 

THD 

t 

Angular rate, Q, (deg/sec) 

FLIGHT PATH ANGLE 

Flight path angle, 7, (deg) 

VEL 

Velocity of nose cap relative to booster (ft/sec) 

XCM 

Aerodynamic moment coefficient 

AA 

Product of reference area times dynamic pressure (lb) 

XM1 

Aerodynamic moment (lb/ft) 

C1P 

Difference between location of CG and center of pressure (ft) 

U1 

Body fixed velocity component in XI direction (ft/sec) 

W1 

Body Fixed velocity component in Z1 direction (ft/sec) 

CN 

Normal force coefficient 

XCG 

X coordinate of CG (ft) 

ZCG 

Z coordinate of CG (ft) 

TH1 

Attitude angle (deg) 

FP 

Force due to pressure difference (lb) 

XQ 

X coordinate of point Q 

ZQ 

Z coordinate of point Q 

XP 

X coordinate of point P 

ZP 

Z coordinate of point P 

WIB 

Velocity component of booster in Z1I direction (ft/sec) 

DELT 

Integration time step (sec) 
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